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Abstract: This work is about pure Yang-Mills theory in four Euclidean dimensions with 
gauge group SU(iV). We study rectangular smeared Wilson loops on the lattice at large 
iV and relatively close to the large- N transition point in their eigenvalue density. We show 
that the string tension can be extracted from these loops but their dependence on shape 
differs from the asymptotic prediction of effective string theory. 
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1. Introduction 

This article is about Wilson loop operators, W(C), in four-dimensional Euclidean SU(iV) 
pure gauge theory. C is a closed, non-self-intersecting continuous curve in M 4 , differentiable 
except at a finite set of points where there are finite discontinuities in the unit tangent 
vector to the curve. W(C) is well defined only after a renormalization that eliminates 
perimeter and corner divergences. 

The inclusion of kinks is a complication one would like to avoid in initial studies. It 
becomes essential if one wants to make the Makeenko-Migdal equations well defined. It 
also is needed if one wishes to employ lattice field theory tools to validate non-perturbative 
assumptions about W(C) for contractible C. 

1.1 Smearing 

A convenient way to renormalize W(C) is to use continuum smearing, henceforth referred to 
as "smearing" . Smearing is a well defined procedure in Euclidean continuum field theory, 
abstracting some more ad-hoc procedures in common use in lattice field theory. It was 
introduced in [1] and a brief review of its history can be found in [2]. Smearing introduces 
an extra parameter, s, of dimension length squared, y/s is an observer's resolution of 
localized objects constructed out of fields. Keeping this resolution nonzero eliminates 
all divergences associated with operator compositeness. Smearing enjoys several useful 
properties: 

• All standard general regularization methods (purely perturbative or lattice) are com- 
patible with smearing. Perturbatively, the counter terms to the classical Lagrangian 
required by ordinary renormalization make all observables constructed out of only 
smeared fields finite - with no restriction on their space-time arguments. Beyond per- 
turbation theory, choosing some reasonable definition of a scale eliminates ultraviolet 
divergences from all smeared observables. 

• 0(4) spacetime invariance remains preserved if the regularization preserves it. 

• Gauge invariance remains preserved if the regularization preserves it. 

• For any C, smearing provides a proper definition of the joint eigenvalue distribution of 
parallel transport round C and these eigenvalues reside on the unit circle. The concept 
of a marginal probability distribution for the parallel transport unitary matrix round 
C in the continuum limit makes sense only after smearing. The classical view of 
parallel transport as an element in the compact SU(./V) group is preserved at the 
quantum renormalized level thanks to smearing, but won't hold with more standard 
methods of continuum regularization. In particular if we allow C to have an exactly 
backtracking segment, subsequent passages will cancel exactly out at the quantum 
level. Smeared Wilson loops obey Polyakov's [3] zigzag symmetry. 
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• At the formal level of the Makeenko-Migdal loop equations, the definition of smearing 
can be extended to loop space [4]. In this formal sense, smearing is defined without 
any reference to the gauge fields and their Lagrangian. The smearing parameter plays 
the role of the evolution variable in a generalization of diffusion to loop space. 

• Small smeared Wilson loops admit a local expansion [2] with well defined "conden- 
sates" at each order. These condensates can be rendered dimensionless by multi- 
plication by a power of s. These dimensionless numbers provide non-perturbative 
definitions of "running" coupling constants at scale s. Unlike their phenomenological 
progenitors [5], [6], the condensates are well defined. They have a calculable per- 
turbative expansion. The terms in this expansion do not vanish. Another way to 
define "condensates" is from the large-momentum structure of a specific observable 
[7]. Such a definition does not ensure that precisely the same condensate enters in 
other observables. The smeared versions are defined in a manner independent of 
the observable. They are dependent on the smearing parameter. This parameter is 
denoted by s in the continuum and by S on the lattice. 

Admittedly, smearing is an artificial device. We do not know of an alternative con- 
struction of a full set of functionals W(C) with properties listed above. 

Preserving the full 0(4) symmetry of space-time in the process of smearing has the 
drawback that smeared correlation functions will not exhibit the unitarity of the underlying 
theory in a transparent manner. For those quantities that have a finite limit as s — > 
transparency will be recovered. One can define a version of smearing which only preserves 
an 0(3) subgroup of 0(4) and thus keep unitarity evident. The cost is the loss of explicit 
Euclidean 0(4) invariance. 

It is important to keep in mind the distinction between the (continuous) smearing we 
are using here and more standard procedures. For example, in [8] smeared rectangular 
Wilson loops were used to extract the interquark potential and ultimately force. The 
smearing was restricted to three directions perpendicular to the "time" direction, which 
was taken to infinity. The smearing was done by iteratively adding one-plaquette-windowed 
paths to the spatial portions of the loop with a weight of 1/2. 

Thus, this older version of smearing differs on two counts from the continuous smearing 
we are employing in this paper: It is only 0(3) invariant and it is defined only on the lattice. 

The proper way to think about this smearing is as a procedure to improve the numerical 
quality of the string energies defined in the asymptotic large-time regime. This smearing 
is a purely lattice technique and plays no role in the subsequent extraction of physically 
relevant parameters. No remnant of this smearing is left in the continuum limit. 

The smearing we use in this work is intentionally designed to make the Wilson loop 
average itself well defined in the continuum, rather than just the interquark force. To 
achieve this an extra, tunable, dimensional parameter is introduced, intuitively describing 
the resolution at which the Wilson loop is observed. 

1.2 Weakly versus strongly coupled regimes 

In 4D SU(iV) pure gauge theory, classical scale invariance is anomalous and gets broken 
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at the quantum level. A scale separating a weakly coupled short distance regime from 
a qualitatively different strongly coupled long distance regime is dynamically generated. 
Observables admit an asymptotic expansion computable in perturbation theory at short 
distances. There are no systematic methods of analytic computation at large distances. 
These two regimes coexist in the same theory and are smoothly connected. The crossover 
is relatively narrow. 

Numerical simulations have established that at large distance the theory confines. This 
is in agreement with experimental data for N = 3. A strong version of the confinement 
postulate is: If we scale a loop C up, C — > pC, W(pC) ~ exp(— p 2 aA m \ n ) as p — > oo with 
A m i n the minimal area of all continuous surfaces in M 4 bounded by C. The string tension a 
has dimension of mass squared and is universal: It does not depend on the shape of C. The 
numerical and the empirical evidence are both restricted to simpler loops. Typically, they 
fit in a plane inside Euclidean M 4 , are non-self-intersecting and not too unusually shaped. 

The string tension can be extracted from unsmeared Wilson loops W(C, s = 0) by 
numerical means. The externally determined resolution scale yfs can be thought of as 
representing an effective thickness of C. For very large loops the fixed finite thickness yfs 
should not enter. One expects W(pC,s) ~ exp(— ap 2 A m \ n {C)) with a independent of s. 

There does not exist so far a mathematically rigorous proof of confinement in the 
continuum limit. Even if we postulate confinement, there is no credible analytical com- 
putation of a in terms of a perturbatively defined scale A. We think methods of effective 
field/string theory could achieve this. In this paper we use the lattice to study in detail the 
crossover from weak to strong coupling. We hope to learn how to quantitatively connect 
the confinement regime to the weakly coupled one and eventually estimate a/A 2 without 
using the lattice in any quantitative way. 

The spirit is the same as in a semiclassical approximation based on instantons [9] . That 
approach was not successful but framed the problem well. It tried to connect the pertur- 
bative regime to a regime in which an MIT bag description held by a crossover described 
using an instanton gas. A more phenomenological approach was based on an extension of 
perturbative OPE to include nonperturbative contributions parametrized by "condensates" 
[5]. This approach was quite successful but is imprecisely defined. It remains unclear how 
a theorist would extract an exact value of a universally meaningful "condensate" even if 
she /he somehow managed to solve QCD exactly in the presence of an acceptable UV cutoff. 

Perhaps it is not by accident that our concrete approach employs an artificial smear- 
ing scale. Smearing provides well defined candidates for SVZ condensates, as already 
mentioned. These condensates also are not directly physical since smearing is somewhat 
ad-hoc. They are universally meaningful though. They would be quantitatively useful only 
if one chose a reasonable level of smearing. A good choice would provide an economical 
parametrization of the short distance - long distance crossover. 

This paper is part of a general strategy. We want to gain analytical control of the 
crossover for simple Wilson loops by exploiting newly established large- N phenomena by 
computer simulation. Then we want to compute a in units of a perturbative scale A 2 . 
Confinement is assumed and one accepts an effective description of the confinement regime 
based on the scale a. This might become practical long before a mathematical proof of 
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confinement is found and without a detailed understanding of what causes it. 

1.3 Large N 

The essence of the weak - strong coupling problem remains present in the limit N — > oo. 
The limit is taken as Ngj^(s) = A(s) is kept fixed [10]. g%{s) is a running coupling constant 
in some standard definition. 

It has been recently established that the weak-strong crossover range collapses into a 
well defined point at infinite N [4]. The transition point depends on the shape of C. As 
the loop C is dilated at fixed non-zero smearing and shape a non-analytic change in the 
single-eigenvalue distribution takes place at a sharply defined scale. For a small loop the 
distribution is insensitive to the compact nature of SU(iV). For a larger loop the full group 
is explored by the parallel transporter round it. For finite iV there is no non-analyticity 
and the full group is felt by parallel transport around all loops. 

Group compactness is a key ingredient for confinement. Perturbation theory is insen- 
sitive to it because it starts from an infinite-range Gaussian integral over YM fields. The 
transition at which the eigenvalues of the parallel transport matrix "discover" the point in 
the group that is farthest from identity is a natural scale for matching perturbation theory 
to a long distance description. Traces of smeared Wilson loops in the fundamental rep- 
resentation remain smooth through the transition even at infinite N, although the single 
eigenvalue distribution is not analytic there. Therefore, traces of smeared Wilson loops 
should match well. The small loop regime is in principle calculable by field theory. More 
recently we have learned that one can also make predictions by analytical means about the 
large loop regime. The framework for doing that is effective string theory. 

1.4 Effective string theory 

Effective string theory [11], [12] bears conceptual similarity to the well known chiral effective 
field theory describing the interactions between soft pions in an SU(iV) gauge theory with 
a moderate number of massless quarks. One assumes that spontaneous chiral symmetry 
breaking occurs via a bilinear condensate and that the finite non-zero pion decay constant 
/tt is the scale typical for effects caused by this breaking. Then, symmetry considerations 
produce a large set of predictions. One has a proof [13] for chiral symmetry breaking at a 
physicist's level of rigor but the proof provides no indication for how to calculate / T in terms 
of A by analytical means. Consider the correlation function of two flavor currents in QCD. 
We have a perturbative description at short distances and a chiral effective field theory 
description at long distances [5], [14]. Joining them at a crossover scale would provide 
some estimates for the ratio f%/A 2 . Resonance contributions come in in the crossover 
regime and the match is complicated. The large-iV limit simplifies matters somewhat 
because the resonances become isolated stable particles coming in as poles. 

The Wilson loop analogy is substantially less developed and we think that time has 
come to look into the problem of matching short to long distances for Wilson loops in some 
detail. The existence of the sharply demarcated matching point on the one hand and the 
smoothness of the observable through this matching point on the other are encouraging. 
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We want to determine by lattice gauge theory methods how an effective string de- 
scription on the strong coupling side of the matching point and close to it works in detail. 
The effective string theory makes predictions for a functional of curves C and the central 
assumption is that these predictions describe Wilson loops. The string tension a is used to 
set the scale in the theory from the outside. The effective string theory in itself does not 
generate any scale. An important issue is how this matching depends on scale- invariant 
features of C. 

The predictions are obtained starting from a limit where the minimal spanning area 
of C is very large. One can ignore any length scale that stays fixed as C is dilated. Very 
large loops are described by the Nambu-Goto action for 2 massless bosonic fields. It is 
impossible to define this limiting theory exactly. This is not needed as an asymptotic 
expansion in inverse loop-size produces well defined terms without a full definition of the 
theory. The predictions made by effective string theory are obtained from the expansion of 
the action around its quadratic approximation. Eventually one reaches an order at which 
non-Nambu-Goto terms are needed. Unlike in the chiral Lagrangian case, this order appears 
to be relatively high. The main ingredient organizing the expansion is the postulated local 
nature of the world-sheet theory. Effective string theory for contractible Wilson loops and 
effective field theory for massless pions differ in scope. The effective chiral Lagrangian 
applies to functions of a finite number of scales, while effective string theory applies to a 
functional of a continuum of scales. 

String theory would require the inclusion of handles in the calculation of corrections. 
It is believed that handles can be neglected at infinite N. Thus, in the 't Hooft limit, 
one ends up using just purely field-theoretical methods of two-dimensional field theory 
when one imposes on the effective string description symmetry restrictions coming from 
the original theory. 

In this numerical work we do not have data of quality needed to identify terms pre- 
dicted by the effective string approach beyond the determinant of Gaussian surface fluc- 
tuations. The leading term states that as the loop is dilated to infinite size, at fixed 
shape, logW(C, s) ~ — o"^4 m i n (C). The two next subleading terms in the asymptotic ex- 
pansion around very large loops come from the determinant of small fluctuations of the 
surface bounded by C about its absolutely minimal area configuration. This configuration 
is assumed unique and well separated from other minima. The first subleading term is pro- 
portional to log( J 4 m i n ) and the second is invariant under scalings of C. Except for simple 
contours, it won't be possible to write down explicit formulas for the second subleading 
term, but, for any specific C, the value of the term can be obtained numerically with relative 
ease. 

The effective string theory cannot make predictions for two terms that are present 
in logW(C,s) and come in between its leading and its subleading predictions. These 
unpredictable terms consist of a perimeter term and a corner term. Both are smearing 
dependent. A potential problem then arises of an "interference" between further smearing- 
dependent subleading terms in W(C, s) and smearing independent terms coming from the 
effective string theory. The consequence of this "interference" is that even given an exact 
formula for W(C,s), coming from SU(iV) field theory, we would not be able to check 
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whether the effective string theory works or not because we would not be able to separate 
out the effective-string prediction. The higher-order effective-string predictions could then 
just "melt away" into the exact expression. 

The field theory produces a log(W(C,s)) which diverges as s — > 0. The divergences 
appear in an asymptotic expansion at fixed C in s as s — > 0. Subtracting the C-dependent 
terms that go as ^ and log K sA 2 would yield a finite expression, ft = 1 at tree level 
in perturbation theory and subsequent terms in the leading log approximation can be 
resummed using the Callan-Symazik equation. This produces a term going like log(log sA 2 ). 

Let us consider the case of a very small smearing parameter first: sa <C 1. It makes 
sense now to just subtract terms that diverge as s — > and then set s = 0. This would 
produce as "pure" a Wilson loop as any other regularization, which does not involve smear- 
ing, would. The power divergence is a perimeter term, well separated from other terms. 
It is clear how to subtract it without introducing finite terms that depend on the curve. 
There are no logarithmic divergences proportional to the perimeter. The subtraction of 
the logarithmic divergence depends on the opening angles of the corners. It is unclear how 
to disentangle the remaining finite parts from the subtraction of the corner divergences 
from unrelated terms dependent on shape features of C. By definition we are considering 
an asymptotic expansion in a parameter describing by how much the Wilson loop has been 
dilated relative to a standard size. The coefficients in the asymptotic expansion in the 
inverse of the dilation parameter are nontrivial functions of the shape of the loop, which 
is kept constant throughout. Effective string theory can be used to provide expressions 
for these coefficients. These expressions consist of functions of loop shape, but not overall 
scale. This particular asymptotic expansion differs from other variants, which also are 
produced by effective string theory. For example, one might extract the interquark force 
from rectangular loops and consider the expansion of this force in the distance between the 
quarks measured in units of the string tension, y/aR. Now the coefficients are just pure 
numbers. 

In our context we are left with an open question as to what effective string theory 
predictions for the dependence on shape parameters of C should be compared to. One 
of our objectives in this paper is to get some guidance on this question from numerical 
simulation. 

2. Rough outline of paper 

We have obtained Monte Carlo estimates for smeared rectangular Wilson loops on a hy- 
percubic lattice at various smearing levels, iV's, couplings, volumes and combinations of 
rectangle sides. The estimates were obtained using a data base of 160 uncorrelated equili- 
brated gauge fields we have distributed on a forty node PC cluster. Each cluster node has 
four cores and a total of 24GB of memory to be able to smear and make measurements on 
four distinct gauge fields simultaneously. 

Wilson loops in all distinct orientations and locations were averaged over for each gauge 
configuration separately. For each set of parameters defining the gauge field action and the 
loop we obtain 160 numbers. The set of these numbers is used for the statistical estimation 
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of various physical parameters. Statistical errors are always determined by jackknife with 
the elimination of one single gauge configuration from the set of 160 at a time. 

We start by extracting the string tension from square loops. The infinite-volume and 
large- A limits on the lattice are dealt with first. Then the string tension is extrapolated to 
its continuum limit. The results and extrapolations are validated against a set of L x L+l 
and L x 2L loops. This analysis is done assuming that the term logarithmic in the area 
is precisely the one predicted by effective string theory but making no assumptions about 
the shape-dependent terms. This is achieved by ensuring that shape-dependent terms play 
no role at this stage. Under the same assumption about the logarithmic area dependence 
a loop-shape dependent number is extracted from the data and compared to the string 
prediction. The data is revisited and using global fits the coefficient of the logarithm of the 
area is determined. We obtain numbers consistent with the effective string theory value. 
This validates the assumption we made earlier when the coefficient was held fixed. The 
dependence on smearing is addressed throughout. 

3. Parameter choices 

We use the standard Wilson single-plaquette action. The standard coupling (3 is incorpo- 
rated into b = 2N2- The N — > oo limit is taken at fixed b. All our gauge configurations 
are on symmetric hypercubes of side V*. V is the total number of sites and the boundary 
conditions are periodic. We wish to use effective string theory to understand the data and 
the structure of the latter is restricted in the continuum by target space 0(4) invariance. 
In order to maintain as much of the latter as possible at the regularized level we work only 
with symmetric volumes. 

For large A there is a bulk transition close to b = 0.360. For substantially smaller 
values of b the system is in a phase disconnected from continuum Yang-Mills theory. We 
also need to maintain b < 0.369 to be sure that spontaneous Z 4 (A) breaking at A = oo [15] 
is avoided on all our volumes, including our smallest, 12 4 . We mainly use the range of 
couplings 0.359 < b < 0.369. These couplings can produce relatively small and fine lattices. 
A judicious exploitation of continuum large- A reduction allows us to always carry out the 
required extrapolation to infinite volume. Continuum large- A reduction is sometimes also 
referred to as partial reduction. It is a conservative version of reduction introduced in [16]. 
Henceforth we shall use the term "reduction" instead of "continuum reduction" . 

We stored statistically independent gauge fields at intervals Ab = 0.001. Satisfac- 
tory statistical independence for our observables is obtained for gauge fields at neighboring 
6's being separated by 500 complete SU(2) updates combined with 500 complete over- 
relaxation passes. 1 The autocorrelation time is equal to about one quarter of this separa- 
tion. The set of A values we use consists of A = 7, 11, 13, 19, 29. The computer time for 
generating a gauge field configuration goes as N 3 V and this is the primary limitation on 
the (A, V) combinations we use. 

A complete SU(2) update consists of sequential updates of ^N(N — 1) SU(2) subgroups. Similarly, a 
complete overrelaxation update consists of a "reflection" of the entire S\J(N) link matrix. 
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Each measurement proceeds after the gauge fields have been smeared. The smearing 
parameter S has mainly been taken in the range 0.2 < S < 0.4. In some cases we have 
data up to S = 0.52. The separation between sequential smearing levels is AS = 0.04. 

The Wilson loops Wn on the lattice are defined by 

W N (L 1 ,L 2 ,b,S,V) = jj(Trl[U l ). (3.1) 

lec 

The product is over the links / in the order they appear when one goes once round C, a 
rectangle of sides Lx )2 . All our fits will be applied to 

w N {L u L 2 , b, S,V) = - log W N (L U L 2 , b, S, V) . (3.2) 

When the loops are square the two variables L\, L 2 are replaced by one L with the under- 
standing that Li ; 2 = L. 



4. Square loops 

4.1 N — > oo and V — > oc limits 

We first want to determine the limit 

lim ( lim w N (L,b,S,V) ) . (4.1) 

N—too yV->oo J 

Numerically this is nontrivial since we need a good level of accuracy on the limit. Fits 
extracting physical parameters are applied to estimates for this limit. 

Once V is larger than some moderate V c (b), large- ./V reduction provides in principle 
a shortcut allowing one to drop the limit V — )■ oo above. At which iV rough convergence 
is attained will depend on V. This requires numerical tests and fits. In addition to N, 
V, and b the magnitude of finite- volume corrections also depends on the observable: the 
larger it is the bigger the finite-volume effects that need to be overcome are. 

We use two different methods to compute the limit (4.1): 

• Method 1) 

At fixed ./V we compute wjv on volumes that are sufficiently large for finite- volume 
effects to be negligible, then we determine Woo(V = oo) by fitting wj^(V = oo) to 

a>i(V = oo) (22(^ = 00) . , , 

w N (V = 00)= Woo (V = 00) + U N2 1 + 2y m ' • (4.2) 

Here, the other arguments of w are omitted for simplicity. All coefficients depend on 
the observable. 

We have evidence (strong for N = 7 and N = 11, not that strong for N = 19 and 
rather weak for N = 29) that volumes V = 24 4 , 18 4 , 14 4 , 12 4 are sufficiently large 
for N = 7, 11, 19, 29, respectively. This statement applies to the specific set of 
couplings and loop sizes we use. The evidence comes from comparing to results on 
other volumes V and/or other values of N or even from trying to extrapolate from 
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lower values of b. The accuracy of our data does not allow to quantify the finite- 
volume effects and we have to settle for something more qualitative. We managed 
to convince ourselves that the finite- volume systematical deviations are smaller than 
our statistical errors. 

• Method 2) 

The second method makes use of large- N reduction. At fixed V, we first take the 
limit iV — > oo of wn(V) by fitting 

■WOO = «-CI0 + + (4-3) 

So long as the center symmetry stays unbroken, there is no volume dependence in 
the infinite- iV theory, i.e., Woo(V) = Woo(V = oo). 

We determine w^V = 12 4 ) from N = 11, 13, 19, 29 [method 2a)] and w^V = 14 4 ) 
from N = 7, 11, 13, 19 [method 2b)]. 

There is little theoretical doubt that reduction indeed holds as a statement about 
N = oo for values of (b, V) in the allowed region. The limit iV — >• oo is unlikely to 
be uniform in (6, V) or in the size of the loop C and in the level of smearing S. If we 
see good fits to a sum of terms decreasing as we know that we have taken into 
account subleading effects that do have a dependence on V. Only then can we trust 
that the leading term is indeed V- independent. 

Figure 1 shows an example for the three different extrapolation methods at b = 0.368, 
S = 0.4, L = 9. Comparing the results we obtain for lirn/v,v-kx) wn{L, b, S, V) in the three 




Figure 1: Plots of w N {L = 9,6 = 0.368,5 = 0.4, V) as a function of 1/iV 2 : V = 12 4 in green, 
V = 14 4 in blue and V = 24 4 (at N = 7), V = 18 4 (at N = 11) in black. Wc display fits for each 
method: 1) by a black solid line, 2a) by a green solid line and 2b) by a blue solid line. Error bars 
are not visible in the plot. 
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different methods 1), 2a), 2b) we obtain credible numbers for Woo(oo). Some comments are 
in order: 

• We obtain reasonable values of xV^dof for the fits (4.2) and (4.3) for 1), 2a), 2b). 
2a) is an exception where we have xV-^dof up to 6 for b < 0.361. We find good 
agreement between the results for limjv,y->oo vjn(V). This agreement is compatible 
with their statistical accuracy of about 0.1%. The worse x 2 /-^dof at 6 < 0.361 likely 
reflect the impact of the N = oo bulk transition in the respective volumes. We keep 
this in mind in subsequent analyses. 

• We cannot determine the coefficients a>2 very accurately. Truncating the expansions 
(4.2) and (4.3) at 0(N~ 2 ) would result in very large x 2 /-^dof, so a 2 cannot be set to 
zero. 

• Including the N = 29 result in the fit (4.3) for V = 12 4 is crucial for 2a) to agree 
with 2b) and 1) for large loops and large b. Including N = 7 in the V = 12 fit 
would require an additional 1/iV 6 correction in (4.3). See Fig. 1 for an example. 
For volumes close to minimal, there is no useful information to be gained about the 
N, V = 00 limit from numbers obtained at low values of N. 

When the lattice size V * is getting close to the critical lattice size L c (b) at which the center 
symmetry brakes, we need to go to higher N's if we want to compute limjv,v->oo wn(L, b, S, V) 
using method 2). The required computation time scales as N 3 V. 2a) is about 1.75 times 
more expensive than 2b) and 1) is about 2.5 times more expensive than 2b). It is hardly 
possible to conclude from 2b) or 2a) alone that the estimates for Woo(V) are reliable. We 
became confident that we have correctly determined lim7v,y->oo wn{V) only after having 
obtained agreeing results from 1), 2a) and 2b). 

4.2 Lattice string tension at infinite N 

At fixed smearing level S and coupling b we use the shorthand notation: 

Woo (L)= lim w N (L,b,S,V). (4.4) 

N,V— >oo 

For square L x L loops, we expect 

Woo (L) + ~ log L 2 = Cl + c 2 L + aL 2 + o(^A . (4.5) 

The log term comes from the determinant of small fluctuations around the minimal area 
configuration in the effective string description. We shall return to it later. For now, its 
presence is just assumed because including it gives good fits while excluding it gives bad 
fits. 

Neglecting corrections of order , we fit 

\ («.(L + 1) - «»(L) + \ log (l + i)) = , (l + 1) + I + O (-^) (4.6) 
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to a straight line as a function of L + g to determine a and C2. See Fig. 2 for some examples. 
Subsequently we use this determination to subtract the area and the perimeter term from 
the data and fit Woo(L) + j logL 2 — ol? — C2L to a constant, c\. We carry out the fits 
in this order because the numerical contribution to wn of aL 2 + C2L substantially exceed 
that of C\. This will be shown later. The numerical contribution of the perimeter term is 
large because it would diverge like s~5 as s — > 0. Allowing too large a perimeter term has 
the negative effect of reducing Wn and consequentially decreasing the signal to noise ratio. 
Reduced smearing increases the statistical error. 

Most 5x5 loops fall into the neighborhood of the large-iV phase transition in the 
eigenvalue spectrum of the Wilson loop matrix for the b and S values we work with. 
Physically smaller loops will have a single-eigenvalue distribution which has a gap around 
-1. We expect confinement to have something to do with the compactness of SU(iV) after 
it is dynamically detected. We know this for a fact in the exactly soluble case of two- 
dimensional YM [17]. In 2D one can analytically separate out all contributions depending 
on eigenvalues exploring the entire unit circle. The remainder of the expression for Wn no 
longer exhibits an area law. Demonstrating the separation requires one to first introduce 
an infrared cutoff because perturbative and non-perturbative contributions mix at infinite 
volume. The absence of confinement has been shown on a sphere, a cylinder and at infinite 
volume by using the Wu-Mandelstam-Leibbrandt IR regularization [18]. 

Therefore, for effective string theory fits we use only loops with L > 6 to determine 
the parameters a, C2, and c\. Wilson loop matrices for such loops have a gapless eigenvalue 
spectrum even in the infinite- N limit. The results presented below are obtained from 
square loops in the range 6 < L < 9. Best fit parameters are obtained by using weighted 
least-square fits. 

Figure 3 shows results obtained for a using the different methods for computing Woo 
for 0.359 < b < 0.369 at smearing level S = 0.4. The three results agree with each other 
within statistical errors. These errors are smallest for method 1). 

Decreasing the smearing level to S = 0.28 results in increasing statistical errors for a. 
Within these errors, o does not depend on S, cf. Tables 1 and 2. This is as we expected 
and indicates that one can use smeared loops as a device to compute the string tension. In 
this sense, the numerical benefit of the type of smearing used before continuous smearing 
was introduced, [8] is maintained. In the older version there is no need to test whether the 
extracted string tension depends on smearing because the infinite time limit projects on 
lowest energy states that by definition are independent of the details of the projector and 
this version of smearing only affects the latter. When one extracts the string tension from 
finite rectangular loops, the independence on the continuous smearing parameter should 
be checked. 

4.3 Continuum limit 

A scale length in lattice units denoted by £ c (6) is used to carry out extrapolations to the 
continuum limit. It is defined at N = 00 using a three- loop calculation of the /3-function 
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for the lattice coupling. The coefficients are written as % = 0, 1, 2. 
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Integrating the RG flow, we define: 



£ c (6) = 0.26 



h 
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Above we have replaced the gauge coupling 6 by 6/(6), 



6/(6)= lim bW N (L 



1, 6,5 = 0, V), 



(4.7) 
(4.8) 

(4.9) 



(4.10) 



(4.11) 



a substitution known as tadpole improvement. In practice this substitution gets one much 
more rapidly into the asymptotic regime where truncated perturbation theory reasonably 
accurately describes the RG flow. The improved infinite- ./V coupling constant 6/(6) is 
defined without smearing. 

The definition of £ c (6) is taken to match with [19]. We only added a numerical pref- 
actor to make £ c (6) ~ L c (b), where L c (b) is given in [15]. This approximation is good to 
10-15% in our range of couplings and would become exact at 6 = oo. It is well known that 
direct continuum extrapolations are subject to large systematic errors since one is too far 
from a truly asymptotic regime. There are many prescriptions for what to use. We do not 



Aw/2 




Figure 2: Plots of 4p = \ {w^L + 1) — w^L) + \ log (l + ^)) obtained with method 1) as a 
function of L + \ at S = 0.4 and b = 0.36 (red), b = 0.362 (black), b = 0.365 (blue) and b = 0.368 
(green). Error bars are not visible in the plot. The straight lines show linear fits through the 
corresponding data points. Only points 6 < L + \ < 9 are used in the fits. The string tension 
values obtained from these fits are collected in Table 1. 
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have a particular opinion about which one is best. Our objective in this paper is not to 
obtain the highest accuracy possible with our data. We preferred to choose a reasonable 
existing prescription to facilitate comparison with other numerical work on closely related 
topics. The net consequence of this is that numerical results for the continuum limit are 
substantially less accurate than numerical results at finite lattice spacing. The work of 
different groups may sometimes be compared only in the continuum limit, and sometimes 
also at finite lattice spacing. It is possible that one ends up with a statistically significant 
discrepancy at finite lattice spacing but no statistically significant discrepancy in the con- 
tinuum limit. This will be the situation for our determination of the string tension, but 
this remark has to be qualified by the comment that one needs to add the observation that 
on the lattice one should get identical numbers for the string tension extracted from the 
two point function of Polyakov loops and that extracted from rectangular loops. 

We separately carry out two two-parameter fits of the relation between the string 
tension a(b) and £ c (6). The two pairs of parameters are denoted by do, d\ and /o, f\: 

and 

^=f^a(b) + f 1 a(b) 2 . (4.13) 

We use ranges 0.359 < b < 0.369 (range A) and 0.362 < b < 0.367 (range B). We also 
use the limited b range (B) since we have observed increasing x 2 's for the infinite- N, V 
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Figure 3: Plots of the string tension values of Table 1: method 1) in black, method 2a) in green, 
method 2b) in blue. 6- values for 2a) and 2b) are slightly shifted to opposite sides of the true b- value 
for visibility. 
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extrapolations using method 2a) for b < 0.361. This was mentioned in Sec. 4.1. Another 
reason is that finite-volume effects increase with increasing b. This reason only applies to 
method 1). 

Results for do and fo using the a(b) values in Table 1 are given in Table 3. Plots of 
the corresponding fit functions (4.12) and (4.13) are shown in Fig. 4. 

The difference between the two fits is a simple indicator of systematic errors induced 
by the truncation of the perturbative series. For all continuum extrapolations do < fo- 
These particular systematic deviations are of the same order as the statistical errors. 

The result of Allton et al. [19] is ->■ 1.85 in our notation with a statistical error of 
1% and a systematic error of 16% on the continuum result. Their numbers were extracted 
from Polyakov loops which are substantially longer than one side of our rectangular loops. 
Apparently, the physical length of string that came into their calculations allows for a 
substantially more accurate determination of the string tension. The systematic errors 
are dominated by the extrapolation to continuum and their relative size is roughly the 
same for us. Although we work at somewhat higher 6-values, this has little impact on the 
extrapolation. Our 6's are still too far from the full set-on of perturbative asymptotics. 

Setting 02 = in the expression for £ c (b) (cf. Eq. (4.10)) results in an increase of about 
26% for do, and an increase of about 29% for fo, cf. Table 4 and Fig. 5. This number is 
too large to commit to a precise estimate of the systematical error. We could be optimistic 
and assume that the next term in perturbation theory would make a smaller correction 
but it is hard to tell. 

Later we shall find some rough estimates of the effective coupling constant at our 
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Figure 4: Plots of a as a function of £~ 2 : method 1) in black, method 2a) in green, method 2b) in 
blue, together with fit functions (4.12) (solid lines) and (4.13) (dashed lines). The fits are obtained 
using 0.359 < b < 0.369 (range A); results for do and fo are given in Table 3. 
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smearing levels and it would come out to be of order unity. The definition of £ c (6) is at 
zero smearing, so these effective couplings are not directly relevant. Nonetheless, the large 
systematic error coming from truncating perturbation theory comes as no surprise. 
In summary, we find that the infinite- A continuum string tension is given by 

lima(&)£ c 2 (6) = 1.6(l)(3). (4.14) 

o— >oo 

The first error is statistical and the second systematic. The systematic error is more of a 
guess than a well founded estimate. Our central number is 2-3 standard deviations smaller 
than that of [19]. In terms of ^-jjg, this translates to: 

a/A^ = 3.4(2)(6). (4.15) 

A previous estimate for the string tension at infinite A extracted from rectangular loops 
has been given in [20]. Expressed in terms of our variables it is cr^c |b/=o.l82 = 2.2(3). The 
results of [20] are claimed to be compatible with [19] at infinite N. These results were 
obtained working at A = 37,47,59 on a 6 4 lattice at small values b = 0.345,0.348,0.350 
assuming large- A" reduction and including folded loops in the analysis. 

While writing up this paper a new study [21] appeared which also deals with rectan- 
gular Wilson loops with the objective to test a new method of full twisted large-A reduc- 
tion [22]. This successful test is carried out on the physical value of the string tension. 
These authors obtain j-^ — = 3.63(3) (statistical error) at A = oo in the continuum if they 

37S 2 

apply the continuum extrapolation method of [19]. This number is fully consistent with 
ours and has very small errors by comparison. The number in [19] is -r-^ — = 3.95(3) (64) 

MS 2 

at A = oo in the continuum. 
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Figure 5: Same as Fig. 4 for /?2 = in £ c (b) (cf. Eq. (4.10)). Fit results arc given in Table 4. 
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There seems to be a disagreement at the statistical level between [19] and our result 
which agrees with that of [21]. The result of [20] seems to side with that of [19], but has too 
large errors to be sure. The systematic error is too big to claim evidence for a difference 
between the string tension extracted from Wilson loops and that extracted from Polyakov 
loop correlators by [19] in the continuum limit. Such a difference would be very difficult 
to accept at the theoretical level. Given the differences between these various simulations 
the case for a real discrepancy between the Wilson loop string tension and the Polyakov 
correlator string tension at the lattice level is not worrying so far. If one ignores smearing, 
these two string tensions ought to be equal already at the lattice level. It might be of 
interest to settle this in future work. 



4.4 String tension at finite N 

We now determine the string tension in lattice units at finite N. Extrapolating to infinite 
iV this would show how the string tension in lattice units converges to its infinite- N limit. 
We do this in order to get a feel about the commutativity of the limits N — > oo and 
UV-cutoff -»• oo. 

We determine the string tension ajy at fixed 6, S, V and N, by fitting C2,n and (Jn in 

2 ("n(L + 1) - w N {L) + \ log (l + £) \ = a N (l + ~) + ^ (4.16) 

to 6 < L < 9 square loop data. We use data obtained on volumes V = 24 4 , 18 4 , 14 4 , 12 4 
for N = 7, 11, 19, 29, respectively, as detailed earlier. We already mentioned in Sec. 4.1 
that for these cases we believe that finite- volume effects are negligible. For N = 7, 11, 19 
we use S = 0.52, and for N = 29 we use S = 0.4. As before, there is no dependence on S 
within statistical errors and the errors decrease with increasing S (cf. Sec. 4.5 below). 

Figure 6 shows plots of the finite- N string tensions as a function of b together with 
the infinite-iV result obtained from Woo by method 1) in Sec. 4.2. For fixed b < 0.367 the 
results for the string tension at finite N are well described by 

a N (b) = a 00 (b) + f ^-. (4.17) 

Results for cr/v and fit parameters (Too, h are given in Table 5. The infinite- N string tension 
obtained in this manner is in good agreement with our previous results in Table 1. Those 
results were determined from Woq. We see that h(b) ~ 10<7oo(&). A relative variation of 
order may seem large. At fixed lattice coupling b one would need to go to N values of 
order 20 to be able to credibly extrapolate to infinite N. 

It will become clear later on, in section 4.6, that the large coefficient of the term is 
replaced by a much smaller one when one looks at a not as a function of b, but as a function 
of b times the plaquette average. The plaquette average also has a finite correction and 
the fact that it works in the manner described is relatively well know, as mentioned in 
Sec. 4.6. 
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4.5 Smearing dependence 



Figures 7 and 8 show plots of wn as a function of L for two different smearing levels 
S = 0.28 and S = 0.52, for N = 11, V = 18 4 , b = 0.365. Also shown are plots of the 
individual terms wjy is composed of (cf. Eq. (5.2)). The variation in ajy is about 1.9%, of 
the order of the statistical errors (1.2% for S = 0.52, 2.6% for S = 0.28). C2 and c\ exhibit 
larger variation. 

These figures allow us to see the relative sizes of the various contributions to the 
Wilson loop at fixed smearing and how these relative contributions change when smearing 
is changed. For a small amount of smearing the negative of the logarithm of the Wilson 
loop is larger (see the red data points and the total fit drawn as a black curve). This 
mainly reflects the larger size of the perimeter term whose contribution is larger than that 
of the area term. The log-term also makes a sizable contribution, smaller than that of 
the area term. The smallest contribution comes from the shape-dependent constant. The 
logarithmic dependence on smearing of the latter has small numerical impact. 

Wilson loops larger than their critical size smoothly merge with their behavior for 
small sizes. One can imagine how all terms except the area one do this. The area term has 
to morph into something else. A likely candidate is a term going as the area squared (for 
planar loops). This term comes from the dimension four condensate which would enter in 
an expansion in loop size small relative to the extent of smearing. For a square loop of side 
I and smearing s, one gets from one gluon exchange wn = ^j^i^fi- So, we only need 
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Figure 6: Plots of the finite- TV string tension 0^(6) from square loops for TV = 7 (blue), TV = 11 
(green), N = 19 (orange), N — 29 (red) together with the infinite- N result (black) obtained by 
method 1), cf. Table 1. 6-valucs for TV = 19 and TV = 29 are slightly shifted in b in the plot for 
visibility. 
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something like l 2 /s < 20 for this approximation to become relevant. 

The level of smearing in lattice units was varied in a large enough range to contain 
levels of smearing in physical scales that are relevant both to the proximity of the large- iV 
transition and to testing whether there is a dependence on smearing in the continuum limit. 
We see that the string tension does not depend on the smearing parameter within statistical 
errors. These errors decrease rapidly with increasing S (cf. Fig. 9 for an example). 

We can check whether the expected continuum divergences as s — > indeed are de- 
tected. On the lattice there will be no divergence as S — > 0, since the lattice spacing 
regulates all divergences, regardless of whether they come from the Lagrangian or from 
the observable. With a reasonable amount of smearing one detects a window where the 
amount of smearing exceeds the lattice spacing influence but is still small enough to exhibit 
the behavior that would have caused a divergence in the continuum. Fig. 10 shows that 
the C2 coefficient of the perimeter term increases linearly with -i, C2 = c^ + /y/S in 
this window. Similarly, Fig. 11 shows that the S-dependence of the L-independent term c\ 
is consistent with a log(5) S — > divergence: c\ = —0.2538 + 0.3278 log S. 

4.6 Continuum limit for finite N 

Most of the relatively large finite- N corrections to the string tension get absorbed when 
considered at fixed improved coupling bi(b,N) rather than fixed b. Figure 12 makes this 
clear. bj(b,N) is determined by the unsmeared plaquette averages at the corresponding 
finite value of TV. Tadpole improvement simultaneously improves the approach to contin- 
uum and to infinite N. Looking at o"jv(6j) instead of (Jn(P) gives a better indication for 
the speed of large- N convergence in the continuum. 
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Figure 7: Plot of w N (L) for N = 11 on V = 18 4 at b = 0.365 and S = 0.28 (red points, error bars 
are not visible), together with analytic functions crL 2 (green, dashed), aL 2 + C2L (red, dotdashed), 
aL 1 + C2L — J logL (blue, dotted), and aL 2 + c-iL — \ logL + c\ (black, solid). The fit parameters 
used to plot the analytic functions were obtained from the data at 6 < L < 9. They are a = 0.02863, 
c 2 = 0.6041, ci = -0.6788. 



- 19 - 



We extrapolate to the continuum at finite N using Eqs. (4.12, 4.13) with the finite- N 
values for £ c , f3 2 and 67. Once we employ the improved coupling scheme, we no longer 
observe a significant ^-dependence in do(N) or fo(N). This is shown in Fig. 13. This 
feature of tadpole improvement [19] was first seen in the context of Polyakov loop correlators 
and is reviewed in [23]. 

4.7 Tree-level continuum perturbation theory 

We calculated the expectation value of a rectangular l\ x l 2 loop in tree-level continuum 
perturbation theory. Diagrammatically, this gives a one-loop integral with a tree-level 
smeared propagator. By "tree-level" we refer to the absence of propagator and vertex 
radiative corrections. 

We parameterize a general closed curve C by x(u) with x(0) = x(l) and x(u) = d u x(u)). 
To leading order in g 2 , wj^- = — log Wjv is given by 



w^{C,s) = \g 2 C 2 ( dux I du 2 D(x(ui) - x(u 2 ); s)5 fiu x fl (ui)x„(u 2 ) , (4.18) 
1 Jo Jo 

D{x]s) = I w^ eW ' x ^ e ~ 2sp2 - (419) 

C 2 is the quadratic Casimir invariant given by C| und = in the fundamental represen- 

tation. We first take a derivative w.r.t. s to bring the integrand into Gaussian form, 

d s D(x-s) = --^^e-^ . (4.20) 

2°7T Z S Z 

Integrating around a rectangular l\ x l 2 loop over u\ >2 then produces error functions. 
Next, we integrate over s, using D(x, s) — > as s — > 00. 

/•oo 

w p N T (h,l 2 ,s) = - ds'd s ,wff(h,l2,s') (4.21) 
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Figure 8: Same as Fig. 7 for S = 0.52. Here, a = 0.02916, c 2 = 0.3813, ci = -0.4624. 
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For s/lf small, we get 
s s 
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(4.24) 



ho(h/h) = ^0(^1/^2) has the following integral representation 
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Beyond the perimeter divergence a loop with kinks will have logarithmic singularities 
as s — > 0. At each kink we denote by 7$ the angle between the tangents. The well known 
expression is: 

9 2 C2 (fundamental) 



N comers 



4vr 2 



-(7,^7,-1)^(^1^). (4.26) 
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Figure 9: String tension as a function of the smearing parameter S for N = 11, b = 0.365, 
V = 18 4 . a is determined using square loops with 6 < L < 9. The horizontal black line corresponds 
to cr(S* = 0.52) = 0.02916. 
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Exactly backtracking segments of C should cancel out for any s / 0. Nonetheless, the 
coefficient of the log(s) divergence above blows up as — > tt. The limits s — > and % — > ir 
do not commute. This is not a surprise, because when a finite loop segment becomes 
exactly backtracking the perimeter changes discontinuously, affecting already the leading 
term in the s — > asymptotic series. 

Note the appearance of the logarithm of the perimeter at leading order in perturbation 
theory. There is a well known fundamental difference between the perimeter divergence 
and the logarithmic one in higher orders of perturbation theory [24] : While the perimeter 
divergence maintains its tree- level dependence on the loop to all orders, the corner diver- 
gence does not because it corresponds to a kink-angle dependent anomalous dimension. 
Were we to choose to define the overall scale as the square root of the minimal area, this 
term would add a shape dependence to the logarithm of the Wilson loop. The logarithm 
of the perimeter comes in from the integral over gluon exchanges where the gluon connects 
point on the opposing sides of the corner. When both endpoints are close to the corner 
the propagator is conformal approximately and this generates a logarithm of the distance 
along each side of the angle. For a very dilated loop having a finite number of separated 
kinks, one expects to leave the conformal regime before any other kink is encountered. 
This would replace the logarithm of the perimeter in the corner divergence term by log t . 
We assumed that s stays fixed, of the order As we mentioned already, there are corner 
terms going as log K s at higher order which sum up to log(logs) in the LLA. 

4.7.1 Perimeter coefficient 

We have seen numerically that there is a term in the logarithm of the smeared Wilson 
loops which is proportional to the perimeter and diverging as s~? for s — > 0. This is a local 
divergence on the loop and therefore ought to be calculable in perturbation theory. One 



Figure 10: Perimeter coefficient C2 (determined from square loops with 6 < L < 9) as a function 
of the smearing parameter S for iV = 11, 6 = 0.365, V = 18 4 . The straight black line shows the fit 
c 2 = -0.2097 + 0.4279/^. 
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expects no log(s) divergent contributions to the perimeter term at all orders in perturbation 
theory. 

From Eqs. (4.22, 4.24) we get the following tree-level formula for the perimeter term 
for a square I x / loop: 

2 (2vr)5 V s 

For example, at b = 0.365 and N = 11, we have obtained c 2 = -0.2097 + 0.4279/v 7 ^ 
(cf. Fig. 10). Matching the coefficient of the L/y/S term with tree- level PT would require 

~ 1.08. This is an indication for the order of magnitude for an effective running 
coupling at this level of smearing. 

4.7.2 Coefficient of log(s) 

The corner divergence at tree level in Eq. (4.24) provides a log(s) term in wj\t: 

g 2 c 2 i 



2 7T 



logs. (4.28) 



For example, at b = 0.365 and N = 11, we obtained 0.3278 for the coefficient of the 
log S term in wj\r (cf. Fig. 11). Matching the numerical result with PT would require 
9-^- ~ 1.03. This is consistent with the perimeter term determination of ^— . 



5. L X 2L loops — shape dependence 

We turn now to a study of the shape dependence of the size-independent term in wm and 
compare it with the effective-string prediction. In general, a shape-dependent parameter 
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Figure 11: L- independent constant c\ (determined using square loops with 6 < L < 9) as a 
function of the smearing parameter S for N = 11, b = 0.365, V = 18 4 . The black line shows the fit 
ci = -0.2538 + 0.3278 log S. 
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characterizing a loop C is a dimensionless number describing C which is invariant under a 
scaling or an M 4 space-time symmetry applied to C. For rectangular L\ x L 2 loops it is 
convenient to introduce the modular invariant shape parameter 



c 



L\ L 2 
L 2 L\ 



(5.1) 



This £ should not be confused with the ^-function that will appear later. 

The accuracy we now need does not permit taking the N — > 00 limit. We restrict our 
attention to the N = 7, 11 data. We shall see that the numbers we compute are identical 
within errors for N = 7 and iV = 11, indicating that it is unlikely that they will change in 
a substantial manner in the N = 00 limit. 

At fixed b, S, V, and fixed finite N, we expect 

1 



w N (L 1 ,L 2 ) + \\ogLiL 2 = ciiv(C) + C2.N Ll t,^ +cf n LiL 2 + O 
1 l 



(5.2) 



Here, the arguments b, S, V are omitted and the single length scale L of squares is replaced 
by Li and L 2 for rectangles. 

We extracted the lattice string tension from square L x L loops at fixed finite N, b, S, 
V in Sec. 4.4. We determined o"jv by fitting the data using Eq. (4.16) with fit parameters 
cat and c 2 ,n- After subtracting area and perimeter terms, we fitted wn(L,L) + \ logL 2 — 
(T^L 2 — c 2: i\[L to a constant. This constant is now denoted by ci i tv(C = 2). 

We now analyze the results obtained for a sequence of rectangular loops at the same 
b, S, V, N with L 2 = 2L\, i.e., C = 3 fixed. Using the results for <jn and c 2 ^ obtained 
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Figure 12: Same data as in Figure 6 now plotted as a function of the improved coupling 6/. The 
solid line shows the result of the continuum extrapolation [method 1) & range A in Table 3] for the 
infmite-iV string tension, a(b) = 1.50/£ c (fe) 2 + 21/£ c (6) 4 , cf. Eq. (4.12). 
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from square loops above, we determine ci.n(( = I) by fitting wn(L,2L) + j log (2L 2 ) — 
0"at2L 2 — C2,at|L to a constant. Figures 14 and 15 show plots of 0x^(2.5) — 01^(2) as a 
function of b for N = 7 and JV = 11. The L-ranges used for fitting ci.tv are 6 < L < 10 for 
square loops and 4 < L < 7 for L x 2L loops. The smallest loop areas included in each set 
are therefore 36 and 32, respectively, putting them close to the large- N phase transition 
point. 

The effective-string prediction for ci(2.5) — ci(2) is 

l log (W))„_ .„ 8664 , ,, 3) 

where rj(x) is the eta-function. We find that the effective-string prediction is smaller than 
the observed values (cf. Figs. 14, 15) by a factor of about 1.5 to 1.7. Within statistical 
errors, our results for 01^(2.5) — ci^(2) do not depend on b, S, or N. 



6. Almost square loops 

We use sequences of almost square loops with sides L\ = L, L<i = L ± 1 to cross check 
our results for the string tension and the shape dependence of c^jy. For these loops, the 
shape-parameter £ changes with L and is given by 

6.1 String tension 

Expanding ci ) at(C) around £ = 2, we obtain from Eq. (4.16) 

X - (w N (L,L + l)-w N (L,L-l) + Ifog^lj =a N L + ^ + ... (6.2) 
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Figure 13: Plot of do(N) (red) and fo(N) (black) determined from finite- N versions of Eqs. (4.12, 
4.13) fitted in the range 0.359 < b < 0.367. 
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We dropped corrections of order jp from the £ expansion and — ^ from the effective-string 
expansion. 

Similarly to the procedure for square loops, we first take lini/vy-^oo Wn and then 
determine the infinite-iV string tension from Woq. Here, we use only method 1), i.e., we 
compute the limit u>oo from data at N = 7, 11, 19, 29 on volumes V = 24 4 , 18 4 , 14 4 , 12 4 , 
respectively (for N = 7 we use V = 24 4 for b > 0.365 and V = 20 4 for b < 0.364). 

To determine a and ci from Eq. (6.2) (at infinite N) we use loops of sizes 6 x 7, 7 x 8, 
8x9. 

The results for the infinite- iV string tension a{b) determined in this manner agree very 
well with those obtained from square loops, cf. Table 6 and Fig. 16. 

6.2 Shape dependence 

Using square and almost square loops, we can study the shape dependence of c\jf{Q. From 
Eq. (5.2), ignoring corrections of order 1/L 4 and 1/crL 4 , we obtain 

1 1 c 'itv( 2 ) 1 
Aw N (L) = w N (L, L) - -w N (L, L + 1) - -w N (L, L - 1) = ^ ^ + • • • (6.3) 

In the first term on the r.h.s c' x N (C) = ■ The second term on the r.h.s. results from 

the |logLiL2 term in Eq. (5.2). To determine c' 1N (2) we multiply Aui^r by L 2 and fit 
to a constant in the range 5 < L < 8 (for an example see Fig. 17). Note the constancy 
of c' t N (2) as a function the gauge coupling. This indicates that the number we extracted 
from the data is already in the continuum limit within quite small errors. 
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Figure 14: Plot of Ci,jv(2.5) - ci,jv(2) for N = 7 (on V 
S = 0.28 (green) S = 0.4 (blue), and S = 0.52 (black). 
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In Awn(L), both the shape-dependent constant ci(C) and the 4log.L1.L2 term result 
in terms of the order 1/L 2 and only their sum can be determined. From the effective-string 
prediction, we would expect L 2 Awn(L) — > 0.0372764 for large L taking into account both 
contributions. This produces the value (^(2) ~ —0.162276. 

The effective string model produces asymptotic predictions for both ci(() and | log L1L2 
from the determinant of Gaussian fluctuations. One may refer to this prediction as a 1- 
loop prediction of effective string theory, to distinguish it from subleading contributions, 
suppressed by powers of the area in units of the string tension. 

Our numerical results for L 2 Awn(L) deviate significantly from the asymptotic string 
prediction (cf. Fig. 17). There is no upfront indication in the data that subleading terms 
in the asymptotic series given by effective string theory play any role since there is no 
dependence on the gauge coupling. Such a dependence would have to appear if the data 
were better described, say, by including a subleading correction. This subleading correction 
would go as one over the area in string tension units and would depend on the gauge cou- 
pling b through the lattice string tension a. Here, we cannot decide whether the deviation 
originates from the shape-dependent constant ci(£) or the logLiL2 term. However, when 
we determine the string tension from square loops, our results seem to be consistent with 
a logL term as predicted by the string model (see also Sec. 7). This indicates that c\(Q is 
responsible for the deviations. 

Taking into account the term coming from the log (i.e., we assume that the co- 
efficient of the logLiL2 term in wn is correctly determined by the string model), our 
results for exceed the asymptotic string prediction by a factor of about 1.6 to 1.8. Since 
ci(2.5) in the string prediction is very well approximated (to an accuracy of 2%) by an 
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Figure 15: Plot of c lfN (2.5) - a, N (2) for N = 11 (on V = 
(red), S = 0.28 (green) S = 0.4 (blue), and S = 0.52 (black). 
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expansion around C, = 2 to linear order, this discrepancy is consistent with the discrepancy 
of ci(2.5) — ci(2) observed in Sec. 5. Our result here is also in agreement within errors 
with [21] who independently report a deviation from string theory. 

We do not observe any significant dependence on b, S or N (cf. Fig. 18 for results at 
S = 0.4). 

There is a question we shall address later on: Is it consistent within effective string 
theory to keep only the leading term? This question is meaningful in the sense that the 
theory predicts a correction of known, universal strength. The expansion in inverse area 
produces a series which likely is divergent. If a subleading correction is larger or equal to 
the leading one, one would conclude that either further terms are needed or, one is outside 
the reach of the asymptotic series in inverse size altogether and no further terms would 
be of any help. Then the fact that the data showed no dependence on the gauge coupling 
b would have to be explained in some other way, unrelated to effective string theory. To 
prepare for this eventuality we need to address two questions first. Do we have numerical 
evidence for the logarithm of area term? What would perturbation theory have to say 
about the Wilson loops in this range? 



6.2.1 Perturbation theory 

Extracting the shape-dependent terms from square I x I and almost square I x (I ±61) loops 
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Figure 16: Plots of the infinitc-JV string tension obtained from rectangular L x L ± 1 loops 
(cf. Table 6) in red, together with results from square loops in black, blue, green (identical to plots 
in Fig. 3). 
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similar to Awn defined in Eq. (6.3) in Sec. 4.7, we obtain 

A^ T (Z, SI, s) = w^ T (l, I, s) - ^w^ T (l, I + SI, s) - ^w^ T (l, I - SI, s) 
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0.260476 and thus 



(6.4) 



A^ T (/, 61, s) = 9 ^ (o.0795774 (jY + o(y 2 ,f 



(6.5) 



We have seen that some terms which diverge at zero smearing, and clearly are outside 
the reach of effective string theory, enter in the logarithm of the smeared Wilson loop in a 
simple additive manner. The case of shape dependence is more complicated and shall be 
discussed later on, when we look at possible explanations for the deviation of the shape 
dependence we measure from the asymptotic prediction of effective string theory. 

7. All loops: validation of the log-term 

So far we have assumed that all our Wilson loops had a prefactor given by (Area) 4. We 
mentioned that attempts to carry out our fits without this term produced substantially 
lower quality fits. We now would like to determine whether the power ^ really is selected 
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Figure 17: Plot of L 2 Aw N (L) for N = 11, V = 18 4 , b = 0.365, S = 0.4 (blue points), result of 
fit to a constant (black solid line, error estimate indicated by the gray band) , and string prediction 
(red dashed line). 
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by our numbers. To do this we need an amount of data and accuracy which does not allow 
us to consider separately different loop shapes or take the infinite- N limit. We do global 
fits to all our data at two values of N. 

At fixed N, b, S, we fit square L x L, almost square L X (L + 1), and rectangular 
L x (2L) loops to 

w N (L 1 ,L 2 ) = oL x L 2 + c 3 log (LiL 2 ) + c 2 Ll I L2 + ci + c 4 (j± + ^ - 2) . (7.1) 

Some results are given in Tables 7 and 8. Going to smaller loops, xV-^dof starts to increase 
significantly. Tables 7 and 8 tell us that the value \ for the exponent of the area is consistent 
with the data. 



7.1 log (Area) versus shape dependence 

We have seen now that one prediction coming from the determinant in the effective string 
description works close to the large- N transition in the eigenvalues and the other does not. 

These two predictions are somewhat different even within effective string theory. The 
determinant of the small fluctuations of the spanning surface around the minimal area one 
is most conveniently evaluated using ^-function regularization. The determinant itself is ill 
defined and ^-function regularization provides one way to extract finite universal features. 
Only such features are conceivably relevant to Wilson loops. 

Within ^-function regularization one has the option to make a decision about how to 
treat the directions perpendicular to the surface which now are fields in the two-dimensional 
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Figure 18: Plot of -c' 1N (2) as a function of b at S = 0.4 for N = 7 on V = 20 4 (red), N = 11 on 
V = 18 4 (green), and N = 19 on V = 14 4 (blue). (We define ci(2) as c[(2) = -L 2 Aw N - |; the 
effective-string prediction is c[(2) w —0.162276.) 
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theory on the world-sheet. This description is supposed to be geometrical without intro- 
ducing any scale. It is convenient to enforce this by thinking in terms of two-dimensional 
gravity. Then, it is natural to view the fields as half-densities [25]. With this conven- 
tion, the power of the area is a regulated number of degrees of freedom, coming from 
C(0; D) where D is the fluctuation operator. The rest of the determinant is just a function 
of the shape parameter £ and comes from the derivative d^(w; D)/dw\ w= o, reflecting the 
eigenvalues of D more directly. 

C(0, D) has additive contributions coming from each kink in our planar curve C: 



y- 



(7-2) 



24tt6» 

Here, the < 9i < 2ir are the angles at each kink measured by an arc contained in the 
interior of C. For a rectangular loop, summing over the two orthogonal directions to the 
surface produces the factor -r. For a backtracking loop the corner term would blow up. For 
a smeared Wilson loop, a backtracking segment makes no contribution. 

Q- function regularization extracts the universal predictions. It is natural to use a two- 
dimensional lattice regularization for the effective string instead. This is so because in the 
strong coupling expansion of the lattice gauge theory one can identify contributions given 
by the exponent of the area of a spanning surface made out of tiles that can be labeled 
by two fields depending on two coordinates on a square world-sheet lattice. It is easy to 
numerically determine the asymptotic expansion in L of the fluctuation determinant for a 
square loop, assuming the two fields to be continuous: 

L-i L-i 

nir mir\ _n 1 



L L J 2 

n=l m=l 



log J [ ( 2 - cos — - cos — J ~ 0.4731L 2 - 0.37645L - - log L - 0.09039 + . . . 

(7.3) 



We see that there is an area term, a perimeter term and a constant but they are absorbed 
into the physical area law, the well defined perimeter and constant terms in the case of 
smeared loops. It is just as easy to do this for a rectangular loop. It is possible to derive the 
asymptotic expansion for rectangular loops by analytical means too. Below we reproduce 
part of Eq. (4.20) from [26]: 

N-1M-1 Ar , 
log ] [ (4 - 2cos ^ - 2cos ^f) ~ MN— - (M + N) log(l + V2) 



n=l m=l 



--logMN+... (7.4) 

Here, G = 1 — 1/3 2 + 1/5 2 ... is Catalan's constant. The universal results using ^-function 
regularization are reproduced. The derivation of [26] shows that the logarithmic term comes 
from modes that vary little (have small wave numbers) in one of the two directions parallel 
to the sides of the rectangle. These modes are the ones most affected by the Dirichlet 
boundary conditions. 
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Looking at the exact expression for the determinant for finite integer (M, N) we find 
that the asymptotic expansion truncated after the constant term provides estimates for the 
logarithm of the determinant at relative accuracy 0.014% for M, N > 5 and 0.0065% for 
M, N > 6. The largest deviation in both sets is at M = N. For loops at fixed aspect ratio 
of 2, M = 2N, the relative accuracy is 0.024% for TV = 4 and 0.0012% for N = 8. In this 
simple case, quite small loops are very well described by the asymptotic series without any 
terms that vanish as M, N — > oo. 

One cannot tell a priori what if anything is left from this numerical observation when 
one looks at real Wilson loops. If this held also for real Wilson loops we would expect 
to see a shape dependence independent on the gauge coupling b for loops as small as we 
looked at. But, then, the number should have agreed with the asymptotic prediction of 
the determinant. 

7.2 Possible explanations of the deviation of the shape dependent constant from 
the prediction of asymptotic effective string theory 

One employs effective string theory under the hypothesis that once the ultimate asymptotic 
regime is entered, it will take complete control of the shape dependence of the functionals 
W(C). Specifically, in the kind of asymptotic expansion in dilatation that we are considering 
(which is different from looking at the separation dependence of the interquark force for 
example) one has, as p — > oo, the following behavior for a dilated loop pC: 



\og(W(pC)) ~ -ap 2 Area min (C) + T P pLength(C) + T^C) log(p) + T 2 (C) 

+ r K (c) + r 3 (c)/ P 2 + r 4 (c)/ P 3 + r 5 (c)/ P 4 + o(i/ P 5 ) . (7.5) 



One can think about p and C as follows: Let the minimal area with some C as boundary 
be unique. Multiply this area by the string tension and get a pure number. Scale C by 
an amount that makes this number equal to 1 and call the so obtained curve C. Choose 
a parametrization of this curve C by a x^(r) with a choice of r such that x 2 = 1. Then, 
the information contained in C is equivalent to the information contained in the set of all 
global 0(4) invariants one can construct out of the function x^. r goes once round C. Note 
that with this convention, the perimeter of pC depends on C. 

C is allowed to have kinks, a is dimensional and has nothing to do with effective 
string theory, except that a > 0. The perimeter coefficient Tp is a non-universal number 
independent of C. Ti(C), ^(C), T^(C), and T$(C) are scale-invariant functions of C and 
universal. T^C) is a scale-invariant function of C with one non-universal overall constant 
multiplicative factor. 



A crucial assumption is that the non-universal function F() is independent of C, and that 
its argument is given by the discontinuity in the tangent of C at the kinks. Without this 
assumption one cannot separate T 2 (C) from Fk(C). With this assumption the Tk term 




(7.6) 



kinks 
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can be eliminated by comparing loops C\ and C2 which have the same set of kinks. In this 
case, effective string theory makes a testable prediction for ^(Ci) — T2(C2)- 

The presence of the terms Tp and Tk in Eq. (7.5) can be motivated in four-dimensional 
Yang-Mills field theory. The Wilson loop has perimeter and corner ultraviolet divergences 
whose removal will introduce some ad-hoc parameters one could not expect effective string 
theory to know about. According to this logic, Wilson loops in three-dimensional Yang- 
Mills theory with continuous gauge groups would not require a Tk term. From the effective 
string theory point of view there is no motivation for such distinction between four and 
three dimensions. 

Equation (7.5) is tested in only a limited manner on rectangular loops. Tk is a fixed 
constant in this set of loops. Since we encountered a deviation at order p° at the level of 
the universality of ^(C), we could suspect that the problem has to do with the presence of 
kinks. Since any kink can be rounded, the requirement on the size of p needed in order to 
justify keeping only the leading terms up to and including T2 in the asymptotic expansion 
in large p involves the local radius of curvature of C, 1Z. For a kink- free C, parametrically 
described by a x(r), we have 



1Z is invariant under re-parameterizations of the boundary and under C — > p C it scales like 
p. With the choice x 2 = 1, we have 1Z = l/|x|, which endows the parameter r with the same 
dimension as that of x. Effectively, in the presence of a single "almost" kink, the Wilson 
loop depends on two hugely disparate scales. There exists a field-theoretical definition of 
the smeared Wilson loop, and we know from experience that observables depending on 
very disparate scales are hard to calculate. It could be though that the effective string 
theory provides a framework which is so different from field theory that this experience is 
irrelevant. 

We now proceed to ask a simpler question: does effective string theory tell us that the 
subleading corrections to the shape dependence are so large that we had no right to compare 
the data to the leading asymptotic prediction? Had we obtained agreement, we probably 
would not have raised this question, like in many previous studies of the interquark force [8], 
which produced agreement in the leading term in the asymptotic expansion relevant to that 
case. 

Within the premise of effective string theory we are working, there is no substantial 
difference for rectangular loops between something as simple as Z2 gauge theory in three 
Euclidean dimensions and SU(oo) pure gauge theory in four Euclidean dimensions. Conse- 
quentially one has ready examples in the recent literature [27] for how to include subleading 
terms. 

In the next subsection we carry out this exercise on our data. 

After that we come back to discuss at a more intuitive level possible differences be- 
tween Z2 three-dimensional gauge theory on the lattice, which is exactly dual to the 
three-dimensional Ising model and hence has a field-theoretical continuum description built 
around the Wilson-Fisher fixed point, and four-dimensional planar SU(iV) gauge theory. 




(7.7) 
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7.3 Subleading terms in effective string theory 

The subleading corrections that have been looked at in detail come in at orders 
corresponding, respectively, to a bulk, a boundary and another bulk correction. Corner 
corrections have not been discussed in the effective string literature, as far as we know. The 
two bulk corrections have universal coefficients, known functions of the shape parameter Q. 
The boundary term has an adjustable coefficient. There exists an unresolved discrepancy 
in the coefficient of p~ 2 : there are two candidates, denoted as C 2 (u) and C2W) differing by 
a u-independent number, where u + u~ l = C,. 2 In the plots we show the contributions of 
each candidate, hoping that one is correct. 3 We have not rechecked the calculations of the 
coefficients by ourselves. We could add these terms as corrections to the effective string 
theory partition function and then take the logarithm or directly to the logarithm. Since 
the numbers differ for our values of loop area, we show both cases. 

We show different forms of contributions up to order p~ 2 in three examples in fig- 
ures 19, 20, and 21. The black points/lines represent the numerical data in all figures. 
With 

w{L u L 2 ) = w(L u L 2 ) + 1 logLxLs - aL x L 2 - c 2 Ll ^ - a(( = 2) (7.8) 

for rectangular L\ x L 2 loops, the colored lines pertaining to the effective string description 
are defined as follows: 

• Red: ^-term coming from the determinant, w(Li, L 2 ) = i log (r)(iu)r](i / u)rj~ 2 (z)) . 

• Solid green: w(Li,L 2 ) = | log (rj(iu)r](i/u)r]~ 2 (i)^ — log ^1 + §j^- % ^ is used to in- 
clude the term of order p~ 2 . 

• Dashed green: w{L\, L 2 ) = | log {r](iu)r](i/u)rj~ 2 (i)) — log ^1 + §j^j^ \ is used to 
include the term of order p~ 2 . 

• Solid orange: w(L\,L 2 ) = ^\og (r]{iu)ri{i / u)r]~ 2 (ij) — ^l^jj is used to include the 
term of order p~ 2 . 

• Dashed orange: w(Li,L 2 ) = g log (r](iu)ri(i/u)r]~ 2 (i)) — §jj^- is used to include the 
term of order p~ 2 . 

For the subleading terms, we use a = 0.02916 (as obtained from square loops for N = 11 
at b = 0.365, cf. Table 5). Note that 

Aw{L) = w(L, L) - 1 w(L, L + 1) - 1 *(L, L - 1) = Aw(L) - 1 log (l - L~ 2 ) . (7.9) 
2 2 8 

2 We use C 2 (u) = (|j) 2 {Au 2 E A {iu) + 2E 2 (iu)E 2 (i/u)) and C 2 (u) = C 2 (u) - ^ as denned in [27]. 

3 We assume a typographical error in the Fourier expansion of the Eisenstein series in [27] (Eq. (A. 4) 
in JHEP 1201, 104) and instead use the expression E 2 k{iu) = 1 + c(i — 2k) ^2n=i ^zk — i {n)e- 2nnu . Here, 
cr m (n) denotes the sum of the m-th powers of the divisors of n. This corresponds to a change by a factor 
of 2 in the definition of u. 
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The adjustable order p~ 3 term is not shown in the figures because we would have to fit its 
coefficient. We also ignore the order p~ 4 term, although it would produce distinguishable 
numbers on the plots 4 

In w(L,2L) —w(L,L) (cf. Fig. 20), square loops are subtracted from rectangular loops 
which have an area twice as large. Therefore, the effective string predictions with C 2 and 
C 2 differ at next-to-leading order. 

It is quite clear that if we wanted to add the missing O (p~ 3 ) boundary term, we could 
get agreement between theory and a few more data points, for smaller loops. One may 
conclude then that effective string theory applies to our data. That the boundary term 
would scale correctly is an automatic consequence of the 6-independence of the data. 

Alternatively, one may simply conclude that for our loops effective string theory makes 
no definite prediction for the shape dependence. We face a standard situation with asymp- 
totic series: adding too many terms is a bad idea and so is having too few. How one looks 
at the data becomes quite subjective. 

Be that as it may, since the main focus of our work is to determine which predictions 
of effective string theory hold close enough to the infinite- N phase transition and therefore 
could come in when one tried to match onto perturbation theory, we are left where we were 
before carrying out this exercise: beyond the area term we can use with some confidence 
also the logarithmic term. Employing effective string prescribed terms of higher order in 
p^ 1 becomes questionable. 

7.4 Rough estimates in perturbative field theory 

We now take an orthogonal view: we try to see how the data could be described by field 
theory in perturbation theory. Just as with the effective string description, we are most 
likely outside the proper reach of this expansion too. It is clear however that one cannot 
dispense with the area term, albeit that it is not predictable by field-theoretical pertur- 
bation theory. Regarding the term logarithmic in the area coming from effective string 
theory, we choose to eliminate it. Perturbation theory will come with its own logarithms 
and there is no objective way to mix them with logarithms coming from effective string 
theory. The relevant lines in the figures 19, 20, 21 are defined below. 

• Solid blue: w(h,l 2 ,s) =w PT (h,l 2 ,s) - (cf. Eqs. (4.22, 4.24)). 

(2tt)2 V V s / 

Note that shifting w by a (Zi /^-independent) constant has no effect on the observables 
we are considering here. 

• Dashed blue: leading term (in the large-Zj/y^ expansion) of the above. 

. Cyan: w(h,l 2 , s) = w^(h,h, s) - (^) + £ log (^)) . 

Replacing the log(s/Z]j2)-term in u; PT (/i, l 2 , s) by log(sA 2 ), this term would no longer 
contribute to the shape dependence. Then the leading order contribution is deter- 
mined by ho(h/h) only (cf. Eq. (4.24)). 

4 We did not manage to reproduce the plots in [27] which include the p~ 4 term (but apparently ignore 
the p -3 term) when we simply implemented the equations therein. We did not pursue this issue because 
we felt that adding more clutter into our figures would be more harmful than informative. 
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Purple: wn = t% T (Zi, k, s 



a 2 c 2 



(2tt)2 



h+k 
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^((OT))- 



This corresponds to the replacement log(s/(Zi + h) ) — > log(sA ) in w (h,h,s). 

For our perturbative estimates shown in Figs. 19, 20, 21, we set the coupling constant to 
^— = 1.03, the value that we obtained from the coefficient of the log(s)-term for b = 0.365, 
N = ll (cf. Sec. 4.7.2). 

Our information from perturbation theory is clearly too limited at this stage to draw 
any concrete conclusions. As far as we went, it seems to work just as well or as badly as 
effective string theory. 

7.5 Suggestions for further research of the shape dependence 

Summarizing the situation so far, it seems that for moderate sizes loops one observes a 
continuum shape dependence which might be explained by perturbation theory and might 
upon extension to much larger loops transit to another value given by effective string theory. 

Taking into account what we know about nonabelian four-dimensional theories, we 
think that the issue of shape dependence deserves further study, albeit somewhat tangential 
to our own long-range project. 

Further numerical checks could be made focusing on the specific issue of shape depen- 
dence. An interesting set of Wilson loops amenable to study on hypercubic lattices have 6{ 
angles equal to The loop is not convex and for such angles the difference between the 
effective-string shape dependence and the field-theoretical one is enhanced. Physically, in 
perturbation theory gluons exchanged between different segments of the loop will predom- 
inantly choose to travel through the "outside" for an obtuse corner angle. On the other 
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Figure 19: Plot of L 2 Aw(L) (black points) for N = 11, V = 18 4 , b = 0.365, S = 0.4 together with 
the result of a fit to a constant (thin dark-grey line, error estimate indicated by the light-gray band). 
The various colored lines show estimates obtained from effective string theory (green, orange) and 
tree-level perturbation theory (blue, purple, cyan) as described in Sees. 7.3 and 7.4. 
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Figure 20: Plot of w(L,2L) - w(L,L) (black points) for N = 11, V = 18 4 , b = 0.365, S = 0.4. 
Numerical results for w are obtained by subtracting area-, perimeter-, and constant-term (with 
coefficients determined from square loop data) from measured w(L\, L2) + \ log(LiI/2). 



hand, surfaces of effective string theory fluctuate above the "inside" of the loop, convex or 
not. In this context it would be also of interest to consider self-intersecting loops given by 
fitting a figure of 8 onto a hypercubic lattice. 

We suggest that the shape dependence of planar Wilson loops presents an interesting 
case for testing the limitations of the effective string approach. We know that a high-energy 




Figure 21: Plot of w(L,2L) (black) and w(L,L) (gray) as a function of the area A for N = 11, 
V = 18 4 , b = 0.365, S — 0.4. Since w(L,L) is zero within errors (except for the smallest loops), 
w(L, 2L) can be viewed as a reasonable estimate for w(\/2A, y/ A/2) — w(VA, VA). The various 
colored lines show effective-string and perturbative predictions for this difference as a function of 
the area A. 
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scattering event in QCD produces after degradation into the IR a pattern of jets that is 
an imprint of perturbation theory. 

• Could it be that even an asymptotically dilated Wilson loop with kinks in Euclidean 
space would have elements of shape dependence that are determined by the field 
theory at short distances and which do not get washed out at large distances? 

• Were that the case, is the effective string theory framework flexible enough to allow 
the inclusion of specific kink terms that can be adjusted to exactly reproduce the 
angle dependence of the anomalous dimensions associated with kinks? 

7.6 How much should one expect 3D lattice Z^ gauge theory to teach us about 
SU(oo) 4D pure Yang-Mills theory? 

In 3D gauge theories with continuous groups there are no corner singularities. The perime- 
ter term is only logarithmically divergent. The renormalization properties of the Wilson 
loop operators are significantly different. We emphasized several features of the corner 
singularities in 4D gauge theory. They have no analogues in 3D. The corner singularities 
play a role in determining the shape dependence of Wilson loops in 4D. Further, the case 
of discrete gauge groups in 3D is substantially distinct from the case of continuous groups. 

One place one can compare the two theories would be on the lattice. For 4D Yang- 
Mills theory one has well defined loop equations. These equations have a formal continuum 
limit in which loops with kinks play a crucial role. The Ising model, which is dual in three 
dimensions to the Z2 gauge model also has a lattice loop equation in terms of a "barbed 
wire" loop boundary [28]. The equation is very different from the four-dimensional lattice 
loop equation for Yang-Mills. 

These equations look stringy, and it has been the folklore that they would lead to an 
exact dual string description. It seems plausible that if such an exact dual description 
exists, the appropriate effective string would bear some closeness to it, say in terms of the 
correct massless degrees of freedom one should use. This may make no difference at leading 
order, but at higher order what is local on the world sheet in one description might be non- 
local in another. The string dual of the Ising model seems not to have a tunable string 
coupling constant, while the string dual of SU(iV) gauge theory seems to have one, which 
can be set to zero by taking N to infinity. There seems to be no analogous freedom and 
limit in the 3D case. Even if handles are exponentially suppressed for the three-dimensional 
Ising case, there is no control on this and no way we know of to actually estimate their 
numerical values; once one works at loop sizes and accuracies sensitive to higher order 
terms in p than the leading one it is difficult to asses how much of a match between data 
and theory one ought to expect. In the large-iV limit of SU(iV) gauge theory there is a 
credible argument that at least one does not have to worry about handle corrections. 

We urge caution in drawing conclusions from 3D Z2 lattice gauge theory about 4D 
Yang-Mills theory. The results of applying effective string theory to the three-dimensional 
Ising model are quite impressive in themselves, without necessarily being relevant in detail 
to four-dimensional gauge theories of the type we have in Nature at a rather fundamental 
level. 
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8. Conclusions 



A stringy parametrization for rectangular Wilson loops holds relatively well all the way 
down to the large-iV transition point. Notably, the scale dependence of the ri(C)log(p)- 
term is consistent with a correct interpretation of the number of degrees of freedom living 
on the worldsheet in the Nambu-Goto case. However, when one gets to the subleading 
level of the dependence on purely shape-dependent parameters, the situation is less clear. 

We are specifically interested in the large-iV limit of four-dimensional SU(iV) pure 
gauge theories and in contours with corner singularities. It could be that there is another 
effective string theory prescription that takes into account the corners in a special way. 
It would describe the large-/? asymptotics of Wilson loops with corners in some different 
manner. To make the Wilson loops we are interested in well defined in the continuum limit, 
one would need to either use smeared loops, or an alternative method, to eliminate the 
ultraviolet divergences inherent in the field-theoretical definition of the observables under 
consideration. Effective string theory is supposed to apply to continuum observables and 
to be free of ultraviolet cutoff effects. After applying constraints resulting from spacetime 
invariance in four dimensions, all ambiguities still present are supposedly parameterizable 
by terms local on the worldsheet, on the boundary and at the corners. These extra terms 
make contributions to the Wilson loops that are suppressed by non-negative powers of the 
inverse area measured in units of the string tension. An effective string theory different 
from the one employed here might use a different set of fields, or allow different kinds of 
additions representing corners, or both. 

So far, it seems possible to try to roughly match rectangular loops across the large- A r 
phase transition. We need to perform more checks on the short scales side of the transition 
to see if one can ultimately turn this into a credible estimate of the ratio cr/A 2 ^. 

Our results might be taken as an indication to consider other observables that admit 
an effective string representation. Much work has been done on the two-point correlation 
function of Polyakov loops. Here, there are no corner divergences to worry about. So far 
one has not established an analogue of the large- ./V phase transition in Wilson loops in this 
case, but we believe this to be possible. So, maybe focusing on Polyakov loop correlations 
would provide a way to temporarily circumvent the issue of dependence on loop shape 
in the presence of kinks. A more esoteric option is to use surface operators where the 
regularization of the operator presents less difficulty [29] . 

In any case, the issue of shape dependence is seen not to be a numerical impediment to 
obtaining a reasonably accurate estimate of aj A 2 by analytical means. Only refinements 
at a later stage might have to take this issue into account. Nonetheless, at a deeper level a 
full understanding of the interplay between field-theoretical properties of kinks and a good 
effective string description of large loops promises to be illuminating. 
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10. Tables 



Method 1) 




Method 2a) 




Method 2b) 


b 


a 


b 


a 


b 


a 


0.359 


0.0492(12) 




0.359 


0.0500(22) 




0.359 


0.0514(22) 


0.360 


0.04597(94) 




0.360 


0.0477(20) 




0.360 


0.0492(21) 


0.361 


0.04170(86) 




0.361 


0.0425(15) 




0.361 


0.0449(16) 


0.362 


0.03946(71) 




0.362 


0.0405(14) 




0.362 


0.0398(15) 


0.363 


0.03767(61) 




0.363 


0.0375(12) 




0.363 


0.0388(13) 


0.364 


0.03544(58) 




0.364 


0.0360(10) 




0.364 


0.0382(10) 


0.365 


0.03266(55) 




0.365 


0.0334(10) 




0.365 


0.0326(11) 


0.366 


0.03105(47) 




0.366 


0.03183(87) 




0.366 


0.03195(98) 


0.367 


0.02878(45) 




0.367 


0.02931(89) 




0.367 


0.02954(77) 


0.368 


0.02701(38) 




0.368 


0.02696(76) 




0.368 


0.02854(76) 


0.369 


0.02508(36) 




0.369 


0.02614(69) 




0.369 


0.02657(72) 



Table 1: Results for the infinite- TV string tension a at S = 0.4. 



Method 1) 




Method 2a) 




Method 2b) 


b 


a 


b 


a 


b 


a 


0.359 


0.0487(22) 




0.359 


0.0480(39) 




0.359 


0.0543(40) 


0.360 


0.0456(17) 




0.360 


0.0486(33) 




0.360 


0.0504(37) 


0.361 


0.0430(16) 




0.361 


0.0434(28) 




0.361 


0.0466(30) 


0.362 


0.0400(13) 




0.362 


0.0434(24) 




0.362 


0.0393(25) 


0.363 


0.0373(11) 




0.363 


0.0363(21) 




0.363 


0.0377(22) 


0.364 


0.03577(99) 




0.364 


0.0347(17) 




0.364 


0.0412(18) 


0.365 


0.03252(84) 




0.365 


0.0331(17) 




0.365 


0.0318(19) 


0.366 


0.03106(73) 




0.366 


0.0314(14) 




0.366 


0.0317(16) 


0.367 


0.02842(77) 




0.367 


0.0288(15) 




0.367 


0.0295(14) 


0.368 


0.02713(62) 




0.368 


0.0271(13) 




0.368 


0.0289(13) 


0.369 


0.02462(57) 




0.369 


0.0255(11) 




0.369 


0.0261(12) 



Table 2: Results for the infinite- N string tension a at S = 0.28. 
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method & range 


d 


X 2 /^dof 


/o 


xV^dof 


1) & A 


1.50(6) 


0.51 


1.62(3) 


0.63 


2a) & A 


1.55(11) 


0.24 


1.66(6) 


0.25 


2b) & A 


1.55(11) 


0.80 


1.68(5) 


0.79 


1) & B 


1.55(14) 


0.42 


1.65(7) 


0.46 


2a) & B 


1.64(27) 


0.14 


1.71(16) 


0.14 


2b) & B 


1.54(26) 


1.56 


1.68(13) 


1.59 



Table 3: Extrapolation to the continuum using a from Table 1. 



fit 


do 


X 2 /^dof 


fo 


X 2 /N doi 


1) & A 


1.89(8) 


0.51 


2.10(3) 


0.69 


2a) & A 


1.94(15) 


0.24 


2.15(7) 


0.26 


2b) & A 


1.95(16) 


0.80 


2.18(7) 


0.79 


1) & B 


1.95(19) 


0.42 


2.13(9) 


0.48 


2a) & B 


2.07(37) 


0.14 


2.20(19) 


0.15 


2b) & B 


1.93(36) 


1.56 


2.17(15) 


1.61 



Table 4: Continuum extrapolations using £ c (6) with setting /?2 = in (4.10). 



b 


CT7 


cm 


C"19 


C29 


Coo 


h 


X 2 /iV d of 


0.359 


0.03413(38) 


0.04337(65) 


0.04726(79) 


0.0480(11) 


0.04927(60) 


-0.741(37) 


0.14 


0.360 


0.03219(33) 


0.03924(57) 


0.04418(71) 


0.04476(87) 


0.04533(53) 


-0.647(33) 


1.3 


0.361 


0.03086(32) 


0.03804(54) 


0.04094(60) 


0.04100(74) 


0.04231(46) 


-0.559(29) 


0.64 


0.362 


0.02904(28) 


0.03546(43) 


0.03839(52) 


0.03888(66) 


0.03974(39) 


-0.524(26) 


0.09 


0.363 


0.02717(27) 


0.03318(40) 


0.03673(44) 


0.03678(54) 


0.03772(34) 


-0.518(23) 


0.91 


0.364 


0.02599(23) 


0.03168(37) 


0.03449(37) 


0.03438(55) 


0.03551(31) 


-0.466(20) 


0.81 


0.365 


0.02420(24) 


0.02916(35) 


0.03124(40) 


0.03240(47) 


0.03259(30) 


-0.412(20) 


0.35 


0.366 


0.02314(21) 


0.02781(30) 


0.03033(38) 


0.03032(41) 


0.03108(27) 


-0.389(18) 


0.68 


0.367 


0.02172(18) 


0.02610(28) 


0.02834(28) 


0.02810(40) 


0.02906(23) 


-0.359(16) 


1.3 


0.368 


0.02032(17) 


0.02468(27) 


0.02677(28) 


0.02647(32) 


0.02740(21) 


-0.345(14) 


2.2 


0.369 


0.01917(16) 


0.02361(25) 


0.02510(25) 


0.02473(31) 


0.02582(20) 


-0.322(13) 


4.6 



Table 5: String tension at finite N = 7, 11, 19, 29, and results of the corresponding extrapolations 
to infinite TV (obtained by fitting and h in Eq. (4.17)). 
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Method 1) for rect. 





(j 


u.ooy 


0.04842(66 J 


0.360 


0.04510(56 J 


U.OOl 


U.U4l0o(0o J 


U.OOz 


U.UOy4l0^4lO J 


0.363 


0.03744(45) 


0.364 


0.03486(39) 


0.365 


0.03239(37) 


0.366 


0.03068(29) 


0.367 


0.02858(29) 


0.368 


0.02688(25) 


0.369 


0.02489(26) 



Table 6: Results for the infinite- N string tension a from rectangular L x L±l loops at S = 0.4. 



b 


S 


a 


C3 


C2 


Cl 


C 4 




0.365 


0.52 


0.02462(29) 


-0.2415(72) 


0.3739(55) 


-0.4423(71) 


-0.2700(40) 


0.25 


0.365 


0.44 


0.02418(33) 


-0.2558(84) 


0.4332(64) 


-0.4895(81) 


-0.2819(46) 


0.31 


0.365 


0.36 


0.02379(38) 


-0.267(10) 


0.5066(75) 


-0.5478(93) 


-0.2918(54) 


0.42 


0.365 


0.28 


0.02348(48) 


-0.273(13) 


0.6028(95) 


-0.622(11) 


-0.2991(67) 


0.45 


0.365 


0.20 


0.02348(73) 


-0.269(19) 


0.734(14) 


-0.714(15) 


-0.3011(100) 


0.77 


0.366 


0.52 


0.02364(29) 


-0.2349(70) 


0.3696(55) 


-0.4414(59) 


-0.2657(40) 


0.20 


0.366 


0.44 


0.02333(33) 


-0.2444(83) 


0.4249(65) 


-0.4906(67) 


-0.2752(47) 


0.15 


0.366 


0.36 


0.02315(39) 


-0.249(10) 


0.4928(78) 


-0.5510(78) 


-0.2816(57) 


0.23 


0.366 


0.28 


0.02307(49) 


-0.249(13) 


0.5832(97) 


-0.6265(98) 


-0.2853(70) 


0.48 


0.366 


0.20 


0.02267(69) 


-0.252(17) 


0.721(13) 


-0.721(14) 


-0.2920(96) 


0.75 


0.367 


0.52 


0.02189(25) 


-0.2377(64) 


0.3775(49) 


-0.4520(60) 


-0.2707(35) 


0.27 


0.367 


0.44 


0.02155(29) 


-0.2488(76) 


0.4331(57) 


-0.4988(68) 


-0.2808(41) 


0.15 


0.367 


0.36 


0.02133(34) 


-0.2565(92) 


0.5018(68) 


-0.5549(80) 


-0.2882(49) 


0.14 


0.367 


0.28 


0.02125(43) 


-0.260(12) 


0.5926(89) 


-0.6226(97) 


-0.2928(63) 


0.20 


0.367 


0.20 


0.02122(64) 


-0.262(17) 


0.724(13) 


-0.706(13) 


-0.2961(91) 


0.24 


0.368 


0.52 


0.02073(24) 


-0.2392(63) 


0.3801(47) 


-0.4529(55) 


-0.2720(34) 


0.34 


0.368 


0.44 


0.02030(27) 


-0.2516(75) 


0.4367(55) 


-0.5004(64) 


-0.2832(40) 


0.19 


0.368 


0.36 


0.01990(33) 


-0.2616(90) 


0.5075(66) 


-0.5582(76) 


-0.2928(48) 


0.16 


0.368 


0.28 


0.01951(41) 


-0.271(11) 


0.6028(84) 


-0.6283(91) 


-0.3017(60) 


0.28 


0.368 


0.20 


0.01905(58) 


-0.280(16) 


0.740(12) 


-0.713(12) 


-0.3107(83) 


0.67 


0.369 


0.52 


0.01996(21) 


-0.2195(57) 


0.3696(42) 


-0.4667(53) 


-0.2637(30) 


0.17 


0.369 


0.44 


0.01968(25) 


-0.2275(70) 


0.4224(50) 


-0.5152(60) 


-0.2726(36) 


0.12 


0.369 


0.36 


0.01945(30) 


-0.2322(86) 


0.4886(62) 


-0.5746(70) 


-0.2794(44) 


0.19 


0.369 


0.28 


0.01921(38) 


-0.236(11) 


0.5793(78) 


-0.6474(82) 


-0.2855(56) 


0.21 


0.369 


0.20 


0.01872(52) 


-0.244(14) 


0.716(11) 


-0.734(11) 


-0.2946(75) 


0.29 



Table 7: Fit results for N = 7 on V — 24 4 . Loop sizes used for the fits are 5x5 to 11x11, 5x6 
to 11 x 12, 4 x 8 to 8 x 16 (which results in iV do f = 14) 
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b 


S 


a 


C3 


C2 


Cl 


c 4 


X 2 /N doi 


0.359 


0.52 


0.04347(77) 


-0.233(15) 


0.334(13) 


-0.421(14) 


-0.2573(95) 


0.41 


0.359 


0.36 


0.04167(99) 


-0.284(20) 


0.505(18) 


-0.540(18) 


-0.294(13) 


0.31 


0.360 


0.52 


0.03912(74) 


-0.262(16) 


0.371(14) 


-0.444(13) 


-0.282(10) 


0.26 


0.360 


0.36 


0.03836(97) 


-0.284(22) 


0.518(18) 


-0.573(19) 


-0.304(13) 


0.28 


0.361 


0.52 


0.03792(75) 


-0.241(16) 


0.355(14) 


-0.442(12) 


-0.2686(100) 


0.08 


0.361 


0.36 


0.03694(98) 


-0.264(22) 


0.503(18) 


-0.573(17) 


-0.293(13) 


0.23 


0.362 


0.52 


0.03531(58) 


-0.260(14) 


0.369(11) 


-0.424(12) 


-0.2756(79) 


0.25 


0.362 


0.36 


0.03456(78) 


-0.281(19) 


0.512(15) 


-0.547(16) 


-0.297(11) 


0.22 


0.362 


0.20 


0.0335(14) 


-0.308(34) 


0.777(27) 


-0.716(26) 


-0.319(19) 


0.16 


0.363 


0.52 


0.03267(55) 


-0.270(13) 


0.384(10) 


-0.435(11) 


-0.2849(75) 


0.41 


0.363 


0.36 


0.03157(73) 


-0.302(18) 


0.533(14) 


-0.552(14) 


-0.312(10) 


0.63 


0.363 


0.20 


0.0319(13) 


-0.289(31) 


0.766(25) 


-0.734(24) 


-0.313(18) 


0.69 


0.364 


0.52 


0.03185(49) 


-0.241(13) 


0.3652(97) 


-0.4440(99) 


-0.2697(70) 


0.12 


0.364 


0.36 


0.03102(66) 


-0.263(17) 


0.505(13) 


-0.564(13) 


-0.2912(94) 


0.17 


0.364 


0.20 


0.0301(13) 


-0.273(33) 


0.758(26) 


-0.750(22) 


-0.309(18) 


0.33 


0.365 


0.52 


0.02929(44) 


-0.257(11) 


0.3815(85) 


-0.4432(96) 


-0.2795(62) 


0.25 


0.365 


0.36 


0.02837(59) 


-0.278(15) 


0.521(12) 


-0.565(13) 


-0.3021(83) 


0.15 


0.365 


0.20 


0.0280(11) 


-0.274(29) 


0.759(22) 


-0.752(20) 


-0.312(16) 


0.31 


0.366 


0.52 


0.02816(46) 


-0.251(11) 


0.3756(87) 


-0.4357(90) 


-0.2731(63) 


0.27 


0.366 


0.36 


0.02745(58) 


-0.276(15) 


0.512(11) 


-0.542(12) 


-0.2940(81) 


0.57 


0.366 


0.20 


0.02689(98) 


-0.287(24) 


0.755(19) 


-0.710(21) 


-0.307(13) 


0.89 


0.367 


0.52 


0.02638(39) 


-0.2388(96) 


0.3763(76) 


-0.4611(85) 


-0.2733(55) 


0.18 


0.367 


0.36 


0.02557(55) 


-0.258(14) 


0.510(11) 


-0.578(11) 


-0.2940(77) 


0.12 


0.367 


0.20 


0.02519(92) 


-0.254(23) 


0.743(18) 


-0.757(19) 


-0.302(13) 


0.31 


0.368 


0.52 


0.02504(37) 


-0.2378(94) 


0.3773(71) 


-0.4614(84) 


-0.2730(52) 


0.12 


0.368 


0.36 


0.02421(50) 


-0.261(13) 


0.5109(98) 


-0.572(11) 


-0.2943(71) 


0.04 


0.368 


0.20 


0.02361(84) 


-0.266(23) 


0.747(17) 


-0.745(19) 


-0.307(12) 


0.38 


0.369 


0.52 


0.02391(34) 


-0.2492(92) 


0.3796(68) 


-0.4335(79) 


-0.2720(49) 


0.29 


0.369 


0.36 


0.02321(45) 


-0.270(12) 


0.5094(90) 


-0.540(11) 


-0.2910(65) 


0.38 


0.369 


0.20 


0.02251(81) 


-0.269(21) 


0.742(16) 


-0.723(17) 


-0.302(11) 


0.45 



Table 8: Fit results for N = 11 on V = 18 4 . Loop sizes used for the fits are 5x5 to 11x11, 5x6 
to 11 x 12, 4 x 8 to 8 x 16. 
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